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Experiences with the multi-level algorithm. 
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Small expectation values are difficult to measure in Monte Carlo calculations as they tend to get swamped 
by noise. Recently an algorithm has been proposed by Liischer and Weisz which allows one to measure 
expectation values which previously could not be measured reliably in Monte Carlo simulations. We will test our 
implementation of this algorithm by looking at Polyakov loop correlators and then explore ways of applying it for 
measuring large Wilson loops. 



1. Introduction 

The physical picture of confinement is often 
given in terms of formation of a flux tube be- 
tween a quark and an anti-quark. This picture is 
also supported by lattice simulations. 

Thinking of the flux tube as a source, the Wil- 
son loop W(r, T) can be considered as correlation 
between sources of length r propagating for time 
T. The theory of this flux tube has been proposed 
as an effective string model 

This is an effective long distance theory and 
therefore to test it, one has to measure physically 
large Wilson loops. The problem however is that 
large Wilson loops have small expectation values. 
That makes it difficult to measure them in simula- 
tions. To circumvent this problem, recent efforts 
to measure properties of the flux tube, have relied 
on asymmetric lattices || . 

Throughout this article, we will work in 3 di- 
mensions and with SU(2). Even in this restricted 
case, (W) ~ f CP 8 will take about 1 year to mea- 
sure to an accuracy of 1% with conventional algo- 
rithms on a lGhz Pill processor. Here wc explore 
a new algorithm which reduces this time to about 
a day. 

2. The Algorithm 

The algorithm is based mainly on the factor- 
ization property of the partition function. This 
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makes it possible to obtain really small expecta- 
tion values as products of expectation values of 
intermediate orders of magnitude. This works for 
the Wilson action to which we restrict ourselves, 
but can be generalized to other actions as long as 
the factorization property holds. We do not go 
into further discussion of the algorithm but just 
present the operational details of our implemen- 
tation of the algorithm. For a detailed discussion 
we refer the reader to Q . 

Suppose we want to measure a Wilson loop of 
extent r and T in x and t directions respectively. 
We choose the time extent of the lattice to be 
a multiple of T and divide the lattice into time 
slices of thickness d. At j3 = 4/c?q = 5 (go is the 
bare coupling), we found d — 2 to be adequate. 
Higher (3 might warrant a higher value of d. 

First we create a 2-link operator (Fig.f A) 

Ta0iS = U* afiUi jS- (f) 

Next take a product of two ^ of these T's to get, 

TctfifS = %, a\i-yvTl nfivS- (2) 

(T)'s are now estimated by an averaging proce- 
dure in which only the time slices are updated 
holding the space-like links at the boundary of 
the slices fixed. They are then multiplied together 
to form the averaged propagator T. Finally the 
sources L\ and L2 are calculated and the full Wil- 
son loop is obtained as (Fig. IB) 

W = L^a^TapysLl p S - (3) 

2 It can be any number depending on T and d 
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Fig.l. Main steps of the multilevel algorithm 
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3. Optimization 

Let us now try to optimize the parameters of 
this algorithm. Probably the most important one 
is the number of sub-lattice updates (iupd) and 
we concentrate on that. Among the others, we 
use three over-relaxation steps to every heat-bath, 
choose the basic time slice as two lattice spacings 
and keep the level of averaging at one. 

For Polyakov loop correlators ((P*P)), the op- 
timal value of iupd was estimated in jjj as the one 
which minimized ytime x (|P*P|). We follow a 
slightly different procedure. The 2-link operator 
T is the main quantity which is affected by the 
sub-lattice averaging. We propose to define the 
norm of this operator as 



N(r) = 




«P(r)W) 2 



(4) 



,a66afi 6aa66a B6a .,„„,« 

Sd "° S5£ --»KS 5 S;5s; 



Fig. 3. error x y/iupd (above) & error 
(below) vs iupd for (W(12, 12)) at /3 = 5. 

size of the loop itself. We estimate the opti- 
mal value of iupd to be the one which minimizes 
N(r) T / 2 x y/iupd. (Note iupd oc time.) N(r) T / 2 
is a number 0((|P*P|)), but measuring N(r) is 
cheaper. We checked that both methods give val- 
ues of iupd which are similar. We also find the 
optimal value of iupd to be reasonably flat and 
not too sensitive to the exact definition of N(r). 

For Wilson loops we have to deal in addition 
with the sources which lie entirely in the space 
slices and whose fluctuation is not reduced by 
this algorithm. Note that it is incorrect to apply 
multi-hit to any of the spatial links in this algo- 
rithm. Fluctuation of the sources greatly reduce 
the optimal number of iupd for Wilson loops. We 
do not have any quantitative estimate on how to 
choose an optimal iupd for the Wilson loop apri- 
ori, but merely point out that if we can estimate 
the order of magnitude of (W) then it is best to 
choose iupd so that each measurement yields a 
value which is of the same order of magnitude. 
In Fig. 3 we plot the error on a 12 x 12 Wilson 
loop against iupd. This plot indicates that the 
optimal value of iupd is ~ 200. In contrast the 
corresponding value for (P*P) is more than 1600. 
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Fig. 2. N(r) versus no. of iupd. 
r= 2 (+), 4 (x), 6 (*) and 8 (□). = 3. 

As seen in Fig. 2, N{r) depends quite strongly 
on r. Moreover the propagator T is a num- 
ber 0{N{r) T l 2 ). Thus the optimal value of 
iupd depends not only on j3, but also on the 



4. Results 

4.1. Polyakov Loop Correlators 

As a test of our program, we compute the string 
tension from (P*P). Neglecting higher eigen- 
states, the static potential between a qq pair is 

V(r)^~]n(P*(r,T)P(0,T)). (5) 
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In the string picture, for large r, V(r) w Vq +ar — 
c/r where Vq is a constant, a is the string tension 
and c (= 7r/24) is a universal number. Next we 
define the force F, and c in terms of V as 
F(r) = [V(r + 1) - V(r - l)]/2 and 
c(r) = [V(r + 1) - 2V(r) + V(r - l)]r 3 /2. 
We tabulate our results below. In each measure- 
Table 1 

V(r), F(r) and c(r) at /3 = 5. Lattice 24 3 . 



r 


V(r) 


F(r) 


c(r) 


2 


0.3460(2) 






3 


0.4634(3) 


0.1129(2) 


-0.1242(8) 


4 


0.5717(5) 


0.1062(2) 


-0.131(2) 


■5 


0.6759(8) 


0.1032(3) 


-0.127(4) 


6 


0.778(1) 


0.1016(4) 


-0.113(9) 


7 


0.879(2) 


0.1009(5) 


-0.08(2) 


8 


0.980(2) 







ment, V(r), F(r) and c(r) were obtained from the 
same sample for all r values to take advantage 
of correlations while computing the differences. 
The averages and errors were obtained using the 
jackknifc method. Finite volume corrections and 
higher state contributions were negligible. 

To get the asymptotic value for a and c, we 
plot F(r) vs 1/r 2 . The slope and intercept of this 
plot gives c and a at r — oo. Wc get c = 0.134(5). 
This is very close to the string value of 7r/24. The 
string tension a^fu = 0.313(3), compares very 
well with Teper's value Q§ of 0.3129(20). 

The above measurement took 10 days on a 1.2 
GHz AMD Athlon processor. 

4.2. Wilson Loops 

Polyakov loops are very good for getting the 
ground state of the flux tube or the string tension, 
but for the excited states and states of higher an- 
gular momentum, it is more useful to evaluate 
Wilson loops j|. 

With the view of computing the spectrum 
of the string, we looked at correlation matrices 
Cij(r,T) between sources ij. As a first step we 
define the source Sj = V (]lfe=i(^fe + a i J2 Sk)), 
with only a single smearing step. Uk is the orig- 
inal link variable, Sk is the sum over all the 
space-like staples and V denotes projection back 
to SU(2). on is the smearing parameter and takes 



values {0.1,0.3,0.5}. 

The energies are extracted by 



E a (r) = —, 



1 



■In 



X a (r,Ti) 



(6) 



T\ - T 2 ' X a (r,T 2 ) ' 

Ti > T2 and A Q (r, T) are the eigenvalues of dj . 
We still do not have enough statistics to get the 
excited states which are highly suppressed, and 
are presently using higher smearing levels and 
better wave functions to extract them. Here we 
present some T dependence for the ground state. 
The apparent downward trend of the data might 

Table 2 

Ground-state for r = 4 determined from various 
values of T x & T 2 . = 5. 



Ti,T 2 


6,4 


12,4 


12,6 


E Q 


0.5728(4) 


0.5723(3) 


0.5722(3) 



be due to the error |6J introduced for using smaller 
values of T\ and T2 ■ The energy determined from 
(P*P), which has T = 24, is 0.5717(5). 

5. Discussion 

With the multilevel algorithm it is possible to 
measure (W) ~ 10~ 8 with 1% error in about a 
day, even on a PC. Even though the algorithm 
works better for Polyakov loops, the improvement 
over multi-hit for (W) < 10~ 6 is quite significant. 
It is possible to couple this algorithm with other 
techniques like smearing to reduce fluctuation of 
the sources and study correlations among them. 
For other applications of the algorithm, see Q]. 
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